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Abstract. The proton capture cross sections of the two most proton rich, stable 
isotopes of Cadmium have been measured for the first time in the energy range relevant 
to the astrophysical p-process between E p = 2.4 and 4.7 MeV. The 106 Cd(p,7) 107 In and 
108 Cd(p,7) 109 In cross sections have been determined using the activation technique. 
Highly enriched as well as natural Cd targets have been irradiated with proton beams 
from both the Van de Graaff and Cyclotron accelerators of the ATOMKI. The cross 
sections have been derived by measuring the 7-radiation following the /3-decay of the 
107 In and 109 In reaction products. The experimental results are presented here and 
are compared with the predictions of Hauser-Feshbach statistical model calculations 
using the NON-SMOKER code. It is found that the calculation reproduces well the 
experimental data. The sensitivity of the model calculations to the proton- and 7- 
strcngths is examined. 
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1. Introduction 

Despite the tremendous experimental and theoretical efforts of recent years, the synthesis 
of the so called p-nuclei (the heavy, proton-rich isotopes which cannot be synthesized 
by neutron capture reactions in the s- or r-process) is still one of the least known 
processes of nucleosynthesis. It is generally accepted that the synthesis of the p-nuclei, 
the astrophysical p-process, involves mainly 7-induced reactions on abundant seed nuclei 
produced at earlier stages of nucleosynthesis by the s- (or to less extent the r-) process 
PQ. During the p-process flow, material from the bottom of the valley of stability is 
driven to the proton rich side by subsequent (7,n) reactions. As the neutron separation 
energy increases while the charged-particle separation energies decrease along this path, 
charged-particle emitting (7, a) and (7,p) reactions become increasingly important for 
more proton-rich nuclei [21 [3j Hj. They can modify significantly the final abundance 
distribution of p-nuclei by feeding isotopic chains with lower charge number Z. 

High energy photons with sufficient energy for dominating photodisintegration 
appear only in explosive nucleosynthetic scenarios. The generally accepted models locate 
the p-process in the deep ONe-rich layers of massive stars either in their pre-supernova 
or supernova phases where temperatures of a few times 10 9 K are reached. 

Comprehensive modeling of the p-process requires on one hand detailed information 
on the stellar environment (temperature, original seed abundances, burning time scale, 
etc.), on the other hand nuclear physics also plays an important role. In p-process 
modeling the reaction rates of hundreds of nuclear reactions involved in nuclear reaction 
networks must be known. All the reaction rates including the dominant 7-induced 
reactions are generally calculated with Hauser-Feshbach type statistical models. The 
rates of 7-induced reactions can be calculated from the inverse capture reactions 
using the detailed balance theorem if the cross sections of the capture reactions are 
known experimentally. While there are compilations for the (n/y) cross sections, very 
few charged particle induced reactions above the iron region have been investigated 
experimentally, leaving the statistical model calculations largely untested. 

Realizing the need for testing experimentally the statistical model calculations in 
the region of proton rich nuclei, several (p,7) and a few (a, 7) reaction cross sections 
have been measured in recent years and the results have been compared with model 
predictions (see e.g. (5] or [6] and references therein, and the KADoNiS database [7]). 
In general, the models are able to reproduce the experimental results within about a 
factor of two, however, some larger deviations are also found. The existing experimental 
database is still not sufficient to check the reliability of model calculations globally, 
therefore further experimental data are highly needed. 

Recently, the (0,7) as well as the (a,n) and (a,p) cross sections of 106 Cd have been 
measured in the energy range relevant to the p-process and compared to statistical model 
calculations |8J. An experiment is also in progress to obtain more information on the 
important a+nucleus optical potential via (a,a) scattering on 106 Cd; preliminary results 
are already available [9]. To give a complete description of 106 Cd from the astrophysical 
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Table 1. Decay parameters of 107 In and 109 In isotopes. Only the strongest 7- 
radiations following the /3-decay of the reaction products which were used for the 
analysis are shown. The data for 107 In are taken from [T2] and for 109 In from [T3] with 
the exception of the 109 In half- life which is taken from a more recent work [13] . 



Product nucleus 


Half life 


Gamma energy [keV] 


Relative intensity 




[hour] 




per decay [%] 


107 In 


0.540 ± 0.005 


204.96 


47.2 ± 0.3 


109 In 


4.168 ± 0.018 


203.50 


73.5 ± 0.5 



p-process point of view, the measurement of its proton capture cross section is also 
necessary. 

In the present work the proton capture cross sections of 106 Cd and 108 Cd (the other 
p-isotope of Cadmium) have been measured. A temperature window of (2.5—3.5) x 10 9 K 
for the astrophysical reaction rates can be covered by our measured energy range. This 
accommodates well the temperatures relevant for the p-process. Consequently, the 
results can be compared with model predictions right at the astrophysically relevant 
energies, no extrapolations are necessary. 

2. Investigated reactions 

Cadmium has 8 stable isotopes with mass numbers 106, 108, 110, 111, 112, 113, 114 
and 116. The two lightest isotopes 106 Cd and 108 Cd are p-isotopes with low natural 
abundances of 1.25% and 0.89%, respectively. Proton-capture on these two isotopes 
leads to unstable In isotopes ( 107 In and 109 In) which decay by /3 + or electron capture 
to 107 Cd and 109 Cd, respectively. For both isotopes the /3-decay is followed by 7- 
radiation which makes it possible to determine the proton-capture cross section using 
the activation method with 7-detection. In this method Cd targets are irradiated by a 
proton beam and the capture cross section is derived from the off-line measurement of 
the decay of reaction products. 

Table 1 shows the decay parameters of the two residual In isotopes. Note that only 
the strongest 7-radiations following the /3-decay of the reaction products are listed. 

Owing to the different decay patterns of the two reaction products, it is possible 
to measure both cross sections in a single activation experiment provided the target 
contains both 106 Cd and 108 Cd isotopes. Cadmium targets of natural isotopic abundance 
could in principle be appropriate for the cross section determination. Proton induced 
reactions on the heavier Cd isotopes, however, can be disturbing if they also lead to 
off-line 7-radiation. In the astrophysically relevant low energy region the cross sections 
of the two investigated capture reactions are very low, thus the elimination of any 
disturbing 7-radiation from the spectra is highly needed. Such disturbing 7-radiation 
can come e.g. from the 110 Cd(p,7) m In or 113 Cd(p,n) 113 In reactions, etc. Therefore, in 
order to avoid the disturbing activity produced by heavier Cd isotopes, enriched targets 
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were used in the lower energy region where the cross section and consequently the yield 
of the studied reactions is very low. 

Both produced In isotopes have metastable states. The proton capture on Cd can 
populate the ground as well as the metastable states in In. In the case of both isotopes, 
however, the metastable states decay exclusively by internal transition to the ground 
state. Moreover, the half-lives of the metastable states are short (54.6 s for 107m In [12] 
and 1.34m for 109m In [13]). Therefore, after a suitable cooling time between the end of 
the irradiation and the beginning of the 7-counting, the metastable states have decayed 
completely to the grounds states and hence the measurement of the ground state decay 
provides information about the total capture cross section. 

3. Experimental procedure 

The experimental setup was similar to the one used for our previous Se(p/y) experiment 
[5], and described in details in a preliminary report on Cd(p/y) [15J. Shortly, thin 
Cd targets have been irradiated with proton beams from both the Van de Graaff and 
Cyclotron accelerators of the ATOMKI. The targets were prepared by evaporating 
natural or highly enriched (enrichment in 106 Cd is 96.47%) metallic Cd onto thin Al 
backing. The thickness of the targets has been measured by weighing. The energy 
range from E p =2.4 to 4.7 MeV was covered in 200 keV steps using the Van de Graaff 
accelerator in the lower and the Cyclotron in the upper part of the energy range. Each 
irradiation lasted about 10 hours with a beam current of 500 nA and the target stability 
was monitored by detecting the backscattered protons from the Cd target. 

Table [2] shows the summary of the irradiations and the list of used targets. Note 
that the targets enriched in 106 Cd contained 2.05% 108 Cd, still higher than the natural 
abundance of 108 Cd (0.89%). 

The induced 7-activity of the targets has been measured with a calibrated HPGe 
detector equipped with 5 cm thick lead shielding. The 7-spectra were taken for at least 
10 hours to follow the decay of both produced In isotopes. The strongest 7-radiations 
from the decay of the two In isotopes are very close to each other in energy (203.5 and 
204.97 keV). However, the energy resolution of the HPGe detector at this low energy is 
about 0.8 keV (FWHM), hence the two peaks could be resolved (see the inset of Fig. [p. 
Moreover, the different half-lives of the two reaction products makes the separation even 
easier. At the beginning of the counting period the 204.96 keV peak from 107 In decay 
dominates, while towards the end of the 10 hour's counting time this isotope has already 
decayed and the 203.5 keV peak form 107 In decay remains. The inset in Fig.[T] shows an 
example where the two peak amplitudes are similar. 

4. Results and conclusion 

The measured cross sections of both investigated reactions cover more than 3 orders of 
magnitude from 3 to 5000 fib. Table [3] lists the experimental results for both reactions 
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Table 2. Details of the irradiations and the targets used. The energy loss has been 
calculated with the SRIM code T6j 



Beam energy 
[keV] 


accelerator 


target type 


target thickness 
[^g/cm 2 ] 


energy loss in target 
[keV] 


O A A A 

2.400 


Van de Graan 


enriched 


595 ± 42 


36.6 


2.600 


Van de Uraan 


enriched 


1 AA _L T 


r a 

5.9 


z.oOO 


Van de Uraan 


enriched 


165 ± 12 


A O 

9.3 


o a c o 
2.952 


cyclotron 


enriched 


43b ± 31 


23.7 


O AAA 
3.000 


Van de Uraan 


natural 


IbO ± 11 


8.1 


O 1 AA 

3.200 


Van de Cjraan 


natural 


17s ± 12 


8. 1 


O 1 AA 

3.200 


Van de Cjraan 


enriched 


1 AA _L T 
1UU ± I 


r o 

5.2 


3.400 


Van de Graaff 


natural 


160 ± 11 


7.5 


3.545 


cyclotron 


enriched 


436 ± 31 


21.2 


3.600 


Van de Graaff 


natural 


178 ± 12 


8.1 


3.744 


cyclotron 


enriched 


146 ± 10 


6.8 


3.942 


cyclotron 


natural 


290 ± 20 


12.5 


4.140 


cyclotron 


natural 


246 ± 17 


10.3 


4.338 


cyclotron 


natural 


246 ± 17 


10.0 


4.536 


cyclotron 


natural 


290 ± 20 


11.4 


4.723 


cyclotron 


natural 


246 ± 17 


9.4 


4.723 


cyclotron 


enriched 


146 ± 10 


5.9 
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Figure 1. Typical activation 7-spectrum taken after the irradiation of a natural Cd 
target with a 3.8 MeV proton beam. This spectrum was taken for 3 hours starting 
8 minutes after the end of the irradiation. The inset shows the two resolved 7-peaks 
from the two reactions studied. The higher energy peaks visible in the spectrum are 
coming from beam-induced activities on heavier Cd isotopes. 
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Table 3. Experimental cross section and S factor of the 106 Cd(p,7) 107 In and 
108 Cd(p,7) 109 In reactions 





106 Cd(p, 7 ) 107 In 






108 Cd(p, 7 ) 109 In 






cross section 


S factor 


F, 

1 - 1 c.m. 


cross section 


S factor 


[MeV] 


[Mb] 


[10 6 MeV-b] 


[MeV] 


[Mb] 


[10 6 MeV-b] 



2, 


,359 


± 





,003 


3.71 ± 


0.43 


201 ± 23 


2.360 


± 


0, 


,003 


3.46 ± 


0.45 


187 


± 


24 


2. 


,573 


± 





,003 


7.96 ± 


0.95 


127 ± 15 


2.573 


± 


0. 


,003 


11.6 ± 


2.0 


186 


± 


32 


2. 


,769 


± 





,003 


21.2 


± 


2.5 


126 ± 15 


2.770 


± 


0. 


,003 


27.3 ± 


3.3 


162 


± 


20 


2. 


,913 


± 


o 


.009 


52.9 


± 


6.1 


162 ± 19 


2.913 




0, 


,009 


70.0 ± 


8.1 


216 


± 


25 


2. 


,968 


± 





,003 


yield too small for analyi 


sis 


2.968 


± 


0, 


,003 


58.0 ± 


7.4 


141 


± 


18 


3 


,166 


± 





,003 


71.5 


± 


9.5 


77.2 ± 10.3 


3.166 


± 


0, 


,003 


157 


± 


19 


170 


± 


21 


3 


,168 


± 





,003 


63.3 


± 


7.3 


67.8 ± 


7.8 


3.168 


± 


0, 


,003 


146 


± 


17 


157 


± 


18 


3 


,364 


± 





,003 


209 


± 


26 


108 ± 13 


3.365 


± 


0, 


,003 


298 


± 


35 


154 


± 


18 


3 


,501 


± 





,011 


223 


± 


26 


72.3 ± 


8.4 


3.505 


± 


0, 


,011 


482 


± 


57 


155 


± 


18 


3 


,562 


± 





,004 


295 


± 


37 


78.4 ± 


9.8 


3.563 


± 


0, 


,004 


391 


± 


45 


104 


± 


12 


3 


,706 


± 





,011 


346 


± 


40 


58.5 ± 


6.8 


3.706 


± 


0. 


,011 


686 


± 


79 


116 


± 


13 


3 


,899 


± 





,012 


528 


± 


64 


50.8 ± 


6.2 


3.900 


± 


0. 


,012 


825 


± 


97 


79.3 


± 


9.3 


4. 


,096 


± 





,012 


936 


± 


114 


52.8 ± 


6.4 


4.097 


± 


0. 


,012 


1862 


± 219 


105 


± 


12 


4. 


,293 


± 





,013 


1232 


± 


149 


42.4 ± 


5.1 


4.293 


± 


0. 


,013 


2744 


± 322 


94.6 


± 


11 


4. 


,488 


± 





,014 


1283 


± 


155 


27.9 ± 


3.4 


4.489 


± 


0. 


,014 


2592 


± 304 


56.5 


± 


6.6 


4. 


,674 


± 





,014 


2329 


± 


280 


33.8 ± 


4.1 


4.675 


± 


0. 


,014 


5392 


± 633 


78.1 


± 


9.2 


4. 


,676 


± 





,014 


1918 


± 


221 


27.7 ± 


3.2 


4.677 


± 


0. 


,014 


4733 


± 545 


68.3 


± 


7.9 



in the form of cross sections as well as astrophysical S-factors. The rows in Table [3] 
correspond directly to the rows of Table [2j There is good agreement between points 
measured at similar energies but with different targets (natural or enriched) or different 
accelerators. 

The quoted uncertainty in the E c m . values correspond to the energy stability of the 
proton beam and to the uncertainty of the energy loss in the target. The uncertainty of 
the cross section values is the quadratic sum of the following partial errors: efficiency 
of the HPGe detector (7%), number of target atoms (6%), current measurement (3%), 
and counting statistics (up to 12%). 

The astrophysical S-factors are plotted in Figs. [2] and [31 They are compared to 
the results of a Hauser-Feshbach statistical model calculation with the NON-SMOKER 
code pm [TT] . It should be noted that the standard calculation is made with global 
parameters which are not adjusted to any local nuclear or reaction properties. There is 
good agreement between experiment and theory across the measured energy range for 
the reaction 106 Cd(p,7) 107 In. There is also good agreement in the case of 108 Cd(p,7) 109 In 
for energies above 3.5 MeV whereas there seems to be a slight systematic underprediction 
below that energy. However, it is still well within the usual 30% error which has 
to be assigned to the global model [T7J [18]. Although it appears as if the data are 
systematically higher below 3.5 MeV in the case of 108 Cd, no correlation with either the 
used targets or the accelerator type was found in a closer inspection of Table El 
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3.5 
E p [MeV] 

Figure 2. Astrophysical S-factor of the 106 Cd(p,7) 107 In reaction as a function of the 
proton cm. energy. The dashed line is the global Hauser-Feshbach statistical model 
prediction of the NON-SMOKER code pi]. 



From Figs. [2] and [3] one might imply that the energy- dependence of the theoretical 
S-factor may be slightly too flat. A stronger increase towards lower energies would 
accommodate the measurements well for both reactions. The calculated cross sections 
and S-factors are mainly dependent on the proton- and 7-strengths. In order to study 
the sensitivity we have performed calculations to study the impact of a variation of 
the strengths by a factor of two. The result is shown in Figs. H] and El The plotted 
sensitivity can assume values between and 1, ranging from no impact to a full factor 
of two change in the cross section and S-factor. It can be seen that the S-factors of both 
reactions are more sensitive to the proton strengths than the 7-strengths across the 
given energy range. This may appear surprising at first glance but it has to be realized 
that the proton widths are smaller than the 7-widths due to the Coulomb suppression 
at such low energies. Consequently, a change in the nuclear properties determining the 
7-strength would only have a limited impact and only at the high end of the measured 
range. The impact would be larger for 108 Cd(p,7) 109 In than for the other reaction. 

The El GDR strength function and the nuclear level density in the compound 
nucleus determine the 7-strength whereas the main property for the proton width is 
the optical potential for protons. The standard calculation shown here uses El and Ml 
transmission coefficients as described in pU] (additionally, E2 transitions were considered 
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E p [MeV] 

Figure 3. Astrophysical S-factor of the 108 Cd(p,7) 109 In reaction as a function of the 
proton cm. energy. The dashed line is the global Hauser-Feshbach statistical model 
prediction of the NON-SMOKER code [II]. 



but they do not contribute significantly here), the nuclear level density of [17], and the 
microscopic proton potential of [19] with low-energy modifications as described in [10J. 
However, since there is reasonable reproduction of the data within the quoted errors it 
is impossible to make a case for a required modification of any of these properties. 

The central quantity important for the application in astrophysical models is the 
astrophysical reaction rate obtained by folding the cross sections with the temperature 
dependent Maxwell-Boltzmann distribution of projectile energies. Because of the limited 
energy range of the data, the relevant integration can only be performed for a limited 
number of temperatures when adopting an upper limit for the error in the integration 
due to the energy cutoff. For the derivation of the astrophysical reaction rates from the 
data, we used the code EXP2RATE [20J which automatically determines the possible 
temperature window and also accounts for the experimental error bars. Allowing for an 
error of at most 15%, we find that the covered temperature range is (2.5 — 3.75) x 10 9 K 
for 106 Cd(p,7) 107 In and (2.5 - 3.5) x 10 9 K for 108 Cd(p,7) 109 In. This covers well the 
relevant p-process temperature range which is usually given as (2.0 — 3.0) x 10 9 K JTJ [3]. 
The resulting rates with their experimental errors are given in Tables H] and [5j 

Comparison of the experimental with the theoretical rates shows that they are 
similar, i.e. compatible with a ratio of unity, within the experimental errors for all given 
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Figure 4. Sensitivity of the 106 Cd(p,7) 107 In cross sections to a variation in the proton- 
and 7-strengths as a function of the proton cm. energy. The sensitivity ranges from 
(no change) to 1 (the cross section is changed by the same factor as the strength). 



Table 4. Astrophysical reaction rates for 106 Cd(p,7) 107 In as function of stellar 
temperature T; only those temperatures are given for which the rates can be computed 
from the data. Also included are the experimental errors. For comparison, also the 
theoretical laboratory rates of [IT] are given, along with the ratio of experimental rate 
to theoretical rate and the stellar enhancement factors SEF. 



T 


exp 


theory 


exp /theory 


SEF 


[10 9 K] 


[cm 3 s -1 mole -1 ] 


[cm 3 s _1 mole -1 ] 






2.50 


1.7797 ± 0.2191 


1.92 


0.93 ± 0.11 


1.01 


2.75 


5.6627 ± 0.7026 








3.00 


(1.5186 ± 0.1893) x 10 1 


1.61 x 10 1 


0.94 ± 0.12 


1.03 


3.25 


(3.5578 ± 0.4445) x 10 1 








3.50 


(7.4729 ± 0.9337) x 10 1 


8.50 x 10 1 


0.88 ± 0.11 


1.05 


3.75 


(1.4346 ± 0.1790) x 10 2 
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Figure 5. Same as Fig. |4]but for the reaction 108 Cd(p,7) 109 In. 
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Table 5. 


Same as Table H but for the reaction 108 Cd(p,7) 109 In. 




T 


exp 


theory 


exp/theory 


SEF 


[10 9 K] 


[cm 3 s _1 mole -1 ] 


[cm 3 s _1 mole -1 ] 






2.50 


2.5625 ± 0.3391 


2.38 


1.08 ± 0.14 


1.00 


2.75 


8.4221 ± 1.1047 








3.00 


(2.3267 ± 0.3040) x 10 1 


2.17 x 10 1 


1.07 ± 0.14 


1.01 


3.25 


(5.5990 ± 0.7313) x 10 1 








3.50 


(1.2044 ± 0.1577) x 10 2 


1.24 x 10 2 


0.97 ± 0.13 


1.00 



temperatures and for both reactions despite the small discrepancies between experiment 
and prediction discussed above. This is due to the fact that the calculation of the 
astrophysical reaction rate involves a weighted average across an appropriate energy 
range and tends to cancel small deviations. 

Finally, a further effect has to be considered under astrophysical conditions. While 
the target nuclei are always in their ground states in a laboratory experiment, nuclei 
in astrophysical plasmas can be thermally excited according to the temperature. This 
has to be taken into account in the calculation of the reaction rate. Experimental and 
theoretical ground state rates are compared in Tables H] and [5j However, the effect of 
thermal excitation is negligible for the cases considered here as is proven by the thermal 
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enhancement factor SEF=r ;(! /r gs which compares the ground state rate r gs to the proper 
stellar rate r* and is also shown in the two tables. 

Summarizing, for the first time we have performed measurements of the two 
reactions 106 Cd(p,7) 107 In and 108 Cd(p,7) 109 In in the energy range relevant for p-process 
nucleosynthesis. S-factors and experimental reaction rates were derived. The results 
excellently confirm the theoretical proton capture rates used in astrophysical models 
for these reactions and thus help to narrow the search for the cause of the remaining 
deficiencies in the nuclear as well as astrophysical p-process modeling. 



Acknowledgments 



This work was supported by OTKA (T42733, T49245, F43408, D48283) and in part 
by the Swiss National Science Foundation (grants 200020-105328, 200020-113984). Zs. 

Fiilop and Gy. Gyiirky are Bolyai fellows. 

[I] M. Arnould and S. Goriely, Phys. Rep. 384 (2003) 1. 

[2] S. E. Wooslcy and W. M. Howard, Astrophys. J. Suppl. 36 (1978) 285. 
[3] T. Rauscher, Phys. Rev. C 73 (2006) 015804. 

[4] W. Rapp, J. Gorres, M. Wiescher, H. Schatz, and F. Kappeler, Asrophys. J. 653 (2006) 474. 

[5] Gy. Gyiirky, Zs. Fiilop, E. Somorjai, M. Kokkoris, S. Galanopoulos, P. Dcmctriou, S. Harissopulos, 

T. Rauscher, and S. Goriely Phys. Rev. C 68 (2003) 055803. 
[6] S. Harissopulos et al, Nucl. Phys. A 758 (2005) 505c. 

[7] I. Dillmann, R. Plag, M. Heil, F. Kappeler, and T. Rauscher, in proceedings of "International 

Symposium on Nuclear Astrophysics - Nuclei in the Cosmos - IX", PoS(NIC-IX)090. 
[8] Gy. Gyiirky et al, Phys. Rev. C 74 (2006) 025805. 

[9] G.G. Kiss, Zs. Fiilop, Gy. Gyiirky, Z. Mate, E. Somorjai, D. Galaviz, A. Kretschmer, K. Sonnabend, 

A. Zilges, Eur. Phys. J. A 27 sOl (2006) 197. 
[10] T. Rauscher and F. K. Thielemann, At. Data Nucl. Data Tables 75 (2000) 1. 

[II] T. Rauscher and F. K. Thielemann, At. Data Nucl. Data Tables 79 (2001) 47. 
[12] J. Blachot, Nucl. Data Sheets 89 (2000) 213. 

[13] D.DeFrenne and E. Jacobs, Nucl. Data Sheets 89 (2000) 481. 

[14] Gy. Gyiirky, Z. Elekes, Zs. Fiilop, G. Kiss, E. Somorjai, A. Palumbo, and M. Wiescher, Phys. Rev. 
C 71 (2005) 057302. 

[15] Gy. Gyiirky G.G. Kiss, Z. Elekes, Zs. Fiilop, E. Somorjai, Eur. Phys. J. A 27 sOl (2006) 141. 

[16] J. F. Ziegler and J. P. Biersack, Code SRIM, Version 2003.20. 

[17] T. Rauscher, F.-K. Thielemann, and K.-L. Kratz, Phys. Rev. C 56 (1997) 1613. 

[18] P. Descouvemont and T. Rauscher, Nucl. Phys. A777 (2006) 137. 

[19] J. P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. C 16 (1977) 80. 

[20; 



T. Rauscher, code EXP2RATE, available at http://download.nucastro.org/codes 



